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oo ■ Abstract: We investigate different metfiods to compute the DC conductance in a quan- j-j-j 

turn wire doped with some impuritied by exploiting the integrability of the theories under 

consideration. As an essential ingredient in all methods we evaluate the reflection and 
> ■ 

. transmission amplitudes of the impurities for a variety of defects. When the impurities in 

0^ , the wire are coupled to an external three dimensional laser field, we predict the generation 

I of harmonic emission spectra. We propose a modified version of the well-known Kubo for- 



CD 



mula, which incorporates the impurities of the system and evaluate the current-current 
' two-point correlation function it involves with the help of a form factor expansion. A 

' comparison with the corresponding quantities computed in a Landauer transport theory 

^ ■ picture is carried out in part II. O 

O 
O 



CJ ■ The work I want to report about is based on a series of papers ^ ^, ^, ^, |5|, |6| with an 

emphasis on the first two. Olalla Castro-Alvaredo will present the second part of this talk. 

^ ■ 1. Generalities on conductance 

In the context of dimensional quantum field theories an impressive arsenal of non- 
perturbative techniques has been developed over the last 25 years. The original motivation 
was to use the lower dimensional set up as a testing ground for general conceptual ideas and 
possibly to apply them in the context of string theory, such that most of the work in this 
area can be characterized very often as rather formal. However, lately the experimental 
techniques have advance to such an extent that one might realistically hope to measure 
various quantities which can be predicted based on these approaches. 

One of those quantities, which is particularly easy to access, is the conductance (con- 
ductivity). It can be measured in general directly without perturbing very much the 
behaviour of the system, e.g. a rigid-lattice bulk metal, such that the uncertainty of ex- 
perimental artefacts is reduced to a minimum. Indeed, there have been some fairly recent 
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measurements [0] of this quantity in 1+1 dimensions and the chahenge is of course to 
explain these data theoretically and possibly inspire more experiments of a similar type. 

There exist two main theoretical descriptions to compute the conductance, the Kubo 
formula pi 0] , which is the outcome of a dynamical linear-response theory and the Landauer- 



Biittinger theory |10|, which is a semi-classical transport theory. The main purpose of the 
work I want to present is a comparison between these two descriptions by employing non- 
perturbative methods of 1+1 dimensional integrable models. It is in this sense the wording 
non-perturbative is to be understood, that is despite the fact that the overall theoretical 
description is of a perturbative nature, within these frameworks we use non-perturbative 
methods. I will concentrate on our proposal of a generalized Kubo formula and in the 
second part, presented by Olalla Castro- Alvaredo, the computations within the Landauer- 
Biittinger transport theory framework will be presented. 

I will start by anticipating the quantities we have to compute. The system we consider 
is a one dimensional quantum wire doped with some impurities (defects). For the time being 
we leave the theory describing the wire and also the nature of the impurities unspecified. 
In linear response theory one essentially needs the Fourier transform of the current-current 
two-point correlation function. This so-called Kubo formula has been adopted to a situation 
with a boundary [^]. Since this only captures effects coming from the constriction of the 
wire a generalization to a set up with defects was needed, which we proposed in [Q] as 

1 r°° 

G"{T) = -]lm^ J_^dte^'^' {J{t)Za ^0))^,^ ■ (1-1) 

Here the defect operator Za enters in-between the two local currents J within the temper- 
ature T and mass m dependent correlation function. The Matsubara frequency is denoted 
by u. 

The other possibility of determining the conductance which we want to study, is a 
generalization of the Landauer-Biittinger transport theory picture. Within this framework 
a proposal for the conductance through a quantum wire with a defect (impurity) has been 
made in |12, 13| 



n- r°° r 

G^{T) = Y, lim 1/ d9 p:{e,T,f^'^\Tt{9)\'-p:{e,T,f^:)\fr{e)\ 



fl.2) 



which we only modify to accommodate parity breaking. This means we allow the trans- 
mission amplitudes for a particle of type i with charge qi passing with rapidity 9 through 
a defect of type a from the left T" (9) and right T" (9) to be different. The density distri- 
bution function pl{9,T,fii) depends on the temperature T, and the potential at the left p,[ 
and right constriction of the wire. 

The main quantities we have to compute before we can evaluate ( |1.1| ) and ( |1.2D are 
the transmission amplitudes Tj, the current-current correlation functions (. . .)T.m and the 
density distributions pi. We obtain all of them non-perturbatively, the T's by means of 



potential scattering theory, e.g. ||lj], the correlation function from a form factor ||l^, 16, 17| 
expansion and the /?'s from a thermodynamic Bethe (TBA) ansatz [18| analysis. 
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2. Impurity systems 

2.1 Constraints from the generalized Yang-Baxter equations 

Let me start with the evahiation of the transmission amphtudes, since they will be required 



in ( |1.1D as well as in {^^). One of the great advantages of integr ability in 1+1 dimensional 
models is that the n-particle scattering matrix factorises into two-particle S-matrices, which 
can be determined by some constraining equations such as the Yang-Baxter |l^ and boot- 
strap equations |2^]. Similar equations hold in the presence of a boundary 22, 23 1 



or a defect [24, ^. It is clear that with regard to the conductance a situation with a 
pure boundary, i.e. non-trivial effects on the constrictions, or purely transmitting defects 
will be rather uninteresting and we would like to consider the case when R and T are 
simultaneously non- vanishing. Unfortunately, it will turn out that for that situation the 
Yang-Baxter equations are so constraining that not many integrable theories will be left 
to consider. Thus this section serves essentially to motivate the study of the free Fermion, 
which after all is very close to a realistic system of electrons propagating in quantum wires. 

We label now particle types by Latin and degrees of freedom of the impurity by Greek 
letters, the bulk scattering matrix by S, and the left/right reflection and transmission am- 
plitudes of the defect by R/R and T/T, respectively. Then the transmission and reflection 
amplitudes are constrained by the "unitarity" relations 

&)R^2{-o) + TL^o)f'^{-e) = (2.1) 



^■',{-o)+n!{e)R)} 

and the crossing-hermiticity relations 



R^^AO)T']{-e)+Til{e)R)U-e) = o, (2.2) 



R'^ie) = R^{-er = Sjji2e)Rf{i7r - 9) , (2.3) 

T^{e) = f^{-ey = f^iin - e) . (2.4) 



The equations ( |2.1D and ( |2.2D also hold after performing a parity transformation, that is 
for <-> i? and T ^ T. 

Depending now on the choice of the initial asymptotic condition one can derive the 
following two non-equivalent sets of generalized Yang-Baxter equations by exploiting the 
associativity of the extended Zamolodchikov-Faddeev algebra |2^, ^ |2^, ^] 

si9i2W® R^iei)]Siei2)[i® R]ii02)] = [i® R^M]Siei2W<^ Rliei)]si9i2), (2.5) 
s{ei2)[i R^M]s{0i2)[i rj(02)] = R^iei) r^(02), (2.6) 

5(0i2)[T,^(e2) Tj(0i)] = [r„^(0i) Tj(02)]5(0i2), (2.7) 



and 



^^(^1) ® R}{02) = RliOi) ® i?^^(02), (2.8) 

[T^{e2) i]s{ei2)[R}{ei) ® i]s{ei2) = rj(02) ® Ri{ei), (2.9) 

[i®f^{e2)]s{ei2)[i®R}{9i)\s{di2) = Ri{di)®r<{e2), (2.10) 

[Tf (0i) 55i]s(^i2)[fj(02) 1] = [i(^f^{e2)]s{ei2)[i®Tj{e{)]. (2.11) 
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We used here the convention {A ® B)^j = A^Bj for the tensor product and abbreviated 
the rapidity sum 612 = 61 + 62 and difference Ou = ^1 — ^2 • Once again the same equations 
also hold for R <^ R and T <^ T. 

Apart from some discrepancies in the indices the equations (^)-( |2l7| ) correspond to a 
more simplified, in the sense that there were no degrees of freedom in the defect and parity 
invariance is assumed, set of equations considered previously in [^]. For diagonal scattering 
it was argued in that one can only have reflection and transmission simultaneously when 
5 = ±1. In Q a more general set up which includes all degrees of freedom was studied. 
A second set of equations ( p.8p -( p!TlD , which is not equivalent to (2.5)-(p77D was found. It 



was shown that in the absence of degrees of freedom in the defect no theory which has 
a non-diagonal bulk scattering matrix admits simultaneous reflection and transmission. 
This result even holds for the completely general case including degrees of freedom in the 
defect upon a mild assumption on the commutativity of R and T in these variables. It 



was further shown that besides S = ±1 also the Federbush model [25| and the generalized 
coupled Federbush models Q allow for i? / and T 7^ 0. 

2.2 Multiple impurity systems 

The most interest situation in impurity systems arises when instead of a single one considers 
multiple defects, since that leads to the occurrence of resonance phenomena and when the 
number of defects tends to inflnity even to band structures. Assuming that the distance 
between the defects is small in comparison to the length of the wire one can easily construct 
the transmission and reflection amplitudes of the multiple defect system from the knowledge 
of the corresponding quantities in the single defect system. For instance for two defects 
one obtains 



T^f^ie) = ^"(f^^'j^^ , RfiO) = Rf{6) + ^^^(^)TnO)Tm (2.12) 

1 - R^\6)Rf{e) 1 - R^{e)Rm 



These expressions allow for a direct intuitive understanding, for instance we note that the 
term [1 - R^ {e)Rf{6)]-^ = E^=i(-Rf (^)^f (6'))" simply resuhs from the inflnite number 
of reflections which we have in-between the two defects. This is of course well known from 
Fabry-Perot type devices of classical and quantum optics. For the case T = T, R = R the 
expressions ( 2.12| ) and ( |2.13|) coincide with the formulae proposed in [^]. When absorbing 



the space dependent phase factor into the defect matrices, the explicit example presented in 
1 24] for the free Fermion perturbed with the energy operator agree almost for T = T,R = R 
with the general formulae ( p. 12 ). They disagree in the sense that the equality of Rf^{6) 



and R'^^{6) does not hold for generic a, (3 as stated in 

It is now straightforward to generalize the expressions for an arbitrary number of 
defects, say n, in a recursive manner 
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Rf{e) 



1 - K 



(0) 



l<k<n. (2.15) 
Similar 



We encoded here the defect degrees of freedom into the vector a={ai,' 
expressions also hold for ff{e) = T;°i-"" (6*) and Rf{d) = Rf'-'^'^ (0). 

Alternatively, we can define, in analogy to standard quantum mechanical methods (see 



e.g. [14|), a transmission matrix which takes the particle i from one side of the defect of 
type a to the other 



-R-{-e)Tt{- 



-Rf{e)T-{e)-^ 



(2.16) 



Then alternatively to the recursive way ( |2.14 ) and ( |2.15| ), we can also compute the multi 
defect transmission and reflection amplitudes as 



^=1 



11 



/ 12 \k=l / 



(2.17) 



^fc=i 



11 



This formulation has the virtue that it is more suitable for numerical computations, since 
it just involves matrix multiplications rather than recurrence operations. In addition it 
allows for an elegant analytical computation of the band structures for n ^ oo, which I 
will however not comment upon further in this talk. 

2.3 Constraints from potential scattering theory 

As we argued in section 2.1., in order to obtain a non-trivial conductance we are lead to 
consider free theories, possibly with some exotic statistics. Trying to be as close as possible 
to some realistic situation, i.e. electrons, we consider first the free Fermion, which with a 



line of defect was first treated in [27|. Thereafter it has also been considered in [g^, 24 1 
and 1 28] from different points of view. In [27, 28, ^ the defect line was taken to be of the 



form of the energy operator and in [29| also a perturbation in form of a single Fermion has 
been considered. In |l[ we treated a much wider class of possible defects. 

Let us consider the Lagrangian density for a complex free Fermion ^ with I defects^ 



\^ 

CD 
O) 

O 
O 



L = ^{i-i^d^ _ m)^ + ^ P"^" (^, V', 5i V^, 5t V)5(x - Xn) 



(2.18) 



n=l 



The defect is described here by the functions P"" (-0, ^, StV') f^tV')) which we assume to 
be linear in the Fermi fields '(/','(/' and their time derivatives. We can now proceed in 



^We use the conventions: 







X 







1 


7 








= {x",x^), p" = (mcosh6»,msmhe), g"" = = e"^ = -e'" = 1 



' ' -10 



5 1 , I i'c, \ J ,10 

1=11, Vc= \ ^(2) ) , tpc= tpLl 



We adopt relativistic units 1 — c = h — niK: e^lSl as mostly used in the particle physics context rather 
than atomic units 1 — e — h — m~ c/137 more natural in atomic physics. 



- 5 - 
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analogy to standard quantum mechanical potential scattering theory (see also |28, p9| ) 
and construct the amplitudes by adequate matching conditions on the field. We consider 
first a single defect at the origin which suffices, since multiple defect amplitudes can be 
constructed from the single defect ones, according to the arguments of the previous section. 
We decompose the fields of the bulk theory as ^p{x) = @{x) Tp-\-{x) + @{—x) ip^{x), with 
@{x) being the Heavyside unit step function, and substitute this ansatz into the equations 
of motion. As a matching condition we read off the factors of the delta function and hence 
obtain the constraints 



27^(V'+(x) - V'-(x))U=o 



dV 



d'ip{x) 



x=0 



d_ 

dt 



dV 



d{dti}{x)) 



(2.19) 



z=0 



We then use for the left (— ) and right (+) parts of ip the well-known Fourier decomposition 
of the free field 



V'f(x) 



with the Weyl spinors 



de 



a,{e)u,{e)e-^^r- + al{e)v,{e)e^^^ 



(2.20) 




-e/2 

,0/2 



(2.21) 



and substitute them into the constraint ( p.l9| ). Treating the equations obtained in this 
manner componentwise, stripping off the integrals, one can bring them thereafter into the 
form 

a,- „(e) = R^{e)aj^^{-e) + r.(0)a,-+(0) , (2.22) 

which defines the refiection and transmission amplitudes in an obvious manner. When 
parity invariance is broken, the corresponding amplitudes from the right to the left do not 
have to be identical and we also have 



(2.23) 



The creation and annihilation operators a|(0) and ai{9) satisfy the usual fermionic anti- 
commutation relations {aj(^i), aj(02)} = 0, {01(6*1), aj(02)} = 27r(5jj(5(^i2). In this way 
one may construct the i?'s and T's for any concrete defect which is of the generic form 
as described in ( |2.18D . After the construction one may convince oneself that the expres- 



sions found this way indeed satisfy the consistency equations like unitarity {2A), ( |2.2D and 

can not be employed for the 



crossing (|2.3D , (p.4|). Unfortunately the equations (|2.l|)-(|2. 
construction, since they are not restrictive enough by themselves to determine the i?'s and 
T's. We consider now some concrete examples: 

2.3.1 Impurities of Luttinger liquid type V^ipjip) = ip{gi +527'^)^' 
Luttinger liquids are of great interest in condensed matter physics, which is one of 
the motivations for our concrete choice of the defect T){'ip,ip) = ip{gi + g2'y^)ip- When 
taking the conformal limit of the defect one obtains an impurity which played a role in this 
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context, see e.g. [31|, after eliminating the bosonic number counting operator. In the way 
outhned above, we compute the related transmission and reflection amplitudes 



Rj{&,9i,92,-y) = Rj{9,gi,g2,y) 

Rji0,gi,g2,-y) = Rj{0,gi,g2,y) 

Tj{9,gi,g2) = T,- (6*, 51, 52) 

Tj{9,gi,g2) = fj{e,gi,g2) 



H92+91 cosh (9) e^'?''"^''^^^ 
{4, + g1 — 52) sinh — 4i{gi + 52 cosh ( 

U{gi -52cosh0)e-2^y™'^''^ii^ 
{4: + gl — gl) sinh 9 — 4:i{gi — g2 cosh i 

(4 + g| -gf)sinhg 

{4 + gf- g'^) sinh 9 - Ai{gi + g2 cosh 9) 

+ -fff)sinh0 

(4 + g(f — 5^2) sinh 9 — 'ii{gi — g2 cosh 9) 



(2.24) 
(2.25) 
(2.26) 
(2.27) 



In the limit lim^j-^o ^(V'; ^) = 5'iV'V') we recover the related results for the T/T's and 
R/R^s for the energy defect operator. For this type of defect we present |Tp and in 
figure 1 with varying parameters in order to illustrate some of the characteristics of these 
functions. 




■ ■ ■ ■ B| = B, = 0.2 y= 10 
- ■B| = b; = 0.2 y = 20 
B, = b; = 0.2 y = 30 



1,5 



1,6 



B| = Bj = 0.2 


y = 


0.25 


- - ■B, = B, = 0.5 


y = 


0.25 


B, = b;= 1.1 


y = 


0.25 





B, = Bj=1.1 y = 0.25 
-B| = ... = 4=1.1 y = 0.25/7 



20 30 40 

E/m = cosh 9 



50 60 20 40 60 80 100 120 140 160 

E/m = cosh 



Figure 1: (a) Single defect with varying coupling constant. jTp and correspond to curves 
starting at and 1 of the same line type, respectively, (b) Double defect with varying distance y . 
(c) Double defect with varying effective coupling constant B = arcsin(— 4gi/(4 + g1)). (d) Double 
defect = dotted line, eight defects = solid line. 
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Part (a) of figure 1 confirms the unitarity relation Part (b) and (c) show the 

typical resonances of a double defect, which become stretched out and pronounced with re- 
spect to the energy when the distance becomes smaller and the coupling constant increases, 
respectively. Part (d) exhibits a general feature, that is when the number of defects is in- 
creased, for fixed distance between the outermost defects, the resonances become more and 
more dense in that region such that one may speak of energy bands. 

2.3.2 The defect V{'il!,i;,dtiJ,dttp) = ig/2{i>dtij - dt^Jip) 

This type of defect reminds on the first non-trivial charge occurring in the free Fermion 
model. In this case we compute by the same means the related transmission and reflection 
amplitudes to 



= R^ie^y) = R^{e,-y) = Rfi9,-y) = _\ , , ^ , (2.28) 



-4ic/cosh6le2*?"""''"^^ 
4ig + tanh 61(4 -|- g'^ cosh^ 6) 



In H] we also computed the T/T's and i?/i?'s for other types of defects, such as 
D = gip^^tp^ T) = gip^^tp, V = g'ip{"/^ ± 'y^)ip ... As an overall conclusion we observed 
that all possible types of parity breaking, that is T 7^ T; R ^ R 01 T ^ T; R = R, etc., 
do occur. We also confirm a general principle one knows well from quantum mechanics, 
namely that parity is preserved when the potential is real, that is in this case the defect 
satisfies V* =1). 

2.4 Impurities coupled to laser fields 

Let us now consider a more complex situation in which a three dimensional laser field hits 
the quantum wire polarized in such a way that it has a vector field component along the 
wire. Since the work of Weyl ||], one knows that matter may be coupled to light by 



means of a local gauge transformation, which reflects itself in the usual minimal coupling 
prescription, i.e. 9^ ^ 5^ — ieA^^, with being the vector gauge potential. The free 
Fermions in the wire are then described by the Lagrangian density 

Ca = i'iil^d^ - m + e-i^'A^)i, . (2.30) 

When the laser field is switched on, we can solve the equation of motion associated to 



- m + e-i^'A^)lP = (2.31) 

by a Gordon- Volkov type solution [^] 



ijjf{x, t) = exp 



/X pt 
dsAi{s,t) {x,t) = ex-p ie / dsAQ{x,s) 'il^j{x,t). (2.32) 



Using now a linearly polarized laser field along the direction of the wire, the vector potential 
can typically be taken in the dipole approximation to be a superposition of monochromatic 
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light with frequency lo, i.e. 

A{t) := Ai{t) = - dsAois) = -- dsE{s) = — / dsf{s)cos{uJs) (2.33) 
^ Jo 2 Jo 2 Jo 

with f{t) being an arbitrary enveloping function equal to zero for t < and t > t, such that 
r denotes the pulse length. In the following we will always take f{t) = 0(t)0(r — t) , with 
Q{x) being again the Heavyside unit step function. The second equality in ( p. 331 ), Aq{x, t) = 



xA{t), follows from the fact that we have to solve (2.32). 

I want to comment on the validity of the dipole approximation in this context. It 
consists usually in neglecting the spatial dependence of the laser field, which is justified 
when xoj < c = 1, where x is a representative scale of the problem considered. In the 
context of atomic physics this is typically the Bohr radius. In the problem investigated 
here, this approximation has to hold over the full spatial range in which the Fermion follows 
the electric field. We can estimate this classically, in which case the maximal amplitude is 
ei?o/u;^ and therefore the following constraint has to hold 

^y = 4t/p<l, (2.34) 

for the dipole approximation to be valid. Due to the fact that a; is a function of lo, we have 
now a lower bound on the frequency rather than an upper one as is more common in the 
context of atomic physics. We have also introduced here the ponderomotive energy Up for 
monochromatic light, that is the average kinetic energy transferred from the laser field to 
the electron in the wire. 

The solutions to the equations of motion of the free system and the one which includes 
the laser field are then related by a factor similar to the gauge transformation from the 
length to the velocity gauge 

ipf{x, t) = exp [ixeA{t)] V'J(x) . (2.35) 

In an analogous fashion one may use the same minimal coupling procedure also to couple 
in addition the laser field to the defect. One has to invoke the equation of motion in order 
to carry this out. For convenience we assume now that the defect is linear in the fields ip 
and tjj. The Lagrangian density for a complex free Fermion ip with i defects T>"{tp,'il^, A^) 
of type a at the position Xn subjected to a laser field then reads 

Cad = Ca + Y, (V^, V-, A^) - x„) . (2.36) 

n=l 

Considering for simplicity first the case of a single defect situated at x = 0, the solution 
to the equation of motion resulting from ( |2.36| ) is taken to be of the form ■ilj^{x,t) = 
Q{x)^pj^_^{x,t) + 0(— x)^j^_(x, i) , which means as before we distinguish here by notation 
the solutions ( |2.35[ ) on the left and right of the defect, V'j^_(x, t) and V'j|^(a^, t), respectively. 
Proceeding as before, the matching condition reads now 



n^(^^(x,t)-e(x,0)U.o=^''""^^'^'^^^ 



dijf{x,t) 



(2.37) 

x=0 
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It is clear, that in this case the transmission and reflection amplitudes will in addition to 
9 and g also depend on the characteristic parameters of the laser field 



T{e,g,Eo,u;,t) and R{9, g, Eo,uj,t) . 



(2.38) 



With regard to the main theme of this talk, it is clear that the laser field can be used to 
control the conductance. For instance defects which have transmission amplitudes of the 
form as the solid line in figure 1 (c), can be used as optically controllable switching devices. 
I want to deviate now slightly from the main line of argument and report briefly on an 
interesting phenomenon one can predict with solutions of the type ( |2.38| ). 

2.5 Harmonic generation 

Let me first briefly explain what harmonics are. The first experimental evidence can be 



traced back to the early sixties [34|. Franken et al found that when hitting a crystalline 
quartz with a weak ultraviolet laser beam of frequency uj, it emits a frequency which is 2uj. 
Generalizing this phenomenon to higher multiples, one says nowadays that high harmonics 
generation is the non-linear response of a medium (a crystal, an atom, a gas, ...) to a laser 
field. Harmonic generation is important, since it allows to convert infrared input radiation 
of frequency u into light in the extreme ultraviolet regime whose frequencies are multiples 
of Lo (even up to order 1000, see e.g. ||3^ for a recent review). A typical experimental 
spectrum is presented in figure 2. 



300000 



In gases, composed of atoms or 
small molecules, this phenomenon is 
well-understood and, to some extent, 
even controllable in the sense that 
the frequency of the highest harmonic, 
the so-called "cut-off" , visible in fig- 
ure 2, can be tuned as well as the in- 
tensities of particular groups of har- 
monics. In more complex systems, 
however, for instance solids, or larger 
molecules, high-harmonic generation 
is still an open problem. This is due 
to the fact that, until a few years 
ago, such systems were expected not 
to survive the strong laser fields one 
needs to produce such effects. How- 
ever, nowadays, with the advent of ultrashort pulses, there exist solid-state materials whose 
damage threshold is beyond the required intensities of lO^^W/cm^ |37]. As a direct con- 
sequence, there is an increasing interest in such materials as potential sources for high- 
harmonics. In fact, several groups are currently investigating this phenomenon in systems 




10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40 

A,(nm) 

Figure 2: Harmonic spectrum for Neon for a Ti:Sa laser 
with A = 795nm. Measured at the Max Born Institut 
Berlin M 



c 

CD 

O 
O 



such as thin crystals ||38|, |39[, carbon nanotubes or organic molecules |4l|, [42|. 



-10- 
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We will therefore try to answer here the question, whether it is possible to generate 
harmonics from solid state devices and as a prototype of such a system we study a quantum 
wire coupled to the laser field in the way described in section 2.4. 

In order to answer that question, we first have to study the spectrum of frequencies 
which is filtered out by the defect while the laser pulse is non-zero. The Fourier transforms 
of the reflection and transmission probabilities provide exactly this information 

T(n,e,Eo,u;,T) = - [ dt\T(e,Eo,u;,t)\'^cos(nt), (2.39) 
^ Jo 
1 f"^ 

n{n,e,Eo,uj,T) = - dt\R{e,Eo,uj,t)\'^cos{nt). (2.40) 
T Jo 

When parity is preserved for the reflection and transmission amplitudes, that is for real 
defects with V* = V, we have |rp + = 1, and it suffices to consider T in the following. 

2.5.1 Type I defects 

Many features can be understood analytically. Taking the laser field in form of monochro- 
matic light in the dipole approximation ( 2.33| ), we may naturally assume that the trans- 



mission probability for some particular defects can be expanded as 

oo 

\Ti{9,Up,u;,t)\^ = Y,t2k{e){^Up)'' sm'"' (Lot). (2.41) 

fc=0 

We shall refer to defects which admit such an expansion as "type I defects" . Assuming that 
the coefficients t2fc(^) become at most 1, we have to restrict our attention to the regime 
4{7p < 1 in order for this expansion to be meaningful for all t. Note that this is no further 
limitation, since it is precisely the same constraint as already encountered for the validity 
of the dipole approximation ( ^.34 ). The functional dependence of ( |2.4l| ) will turn out to 



hold for various explicit defects considered below. Based on this equation, we compute for 
such type of defect 

S TSjnt,P-(!!/2c.)^l 

It is clear from this expression that type I defects will preferably let even multiples of the 
basic frequency uj pass, whose amplitudes will depend on the coefficients t2fc(^)- When we 
choose the pulse length to be integer cycles, i.e. r = 2Trn/oj for n € Z, the expression in 
(2.42) reduces even further. The values at even multiples of the basic frequency are simply 



rji2nu;,e,Up) = (-1)- £ t2fc(^) (f/p)' , _ ' (2-43) 



fe=o 




which becomes independent of the pulse length r. Notice also that the dependence on 
Eq and uj occurs in the combination of the ponderomotive energy Up. Further statements 
require the precise form of the coefficients t2k{0) and can only be made with regard to a 
more concrete form of the defect. 
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2.5.2 Type II defects 

Clearly, not all defects are of the form ( 2.41| ) and we have to consider also expansions of 
the type 



^ Tp2k+p 

\Tn{e,E,/e,u:,tt = ^ cosP(a;t) sin^'^M. (2.44) 

/c,p=0 

We shall refer to defects which admit such an expansion as "type II defects" . In this case 
we obtain 



oo p 



Tn{n,e,Eo/e,u;,T) = 

k,p=0 1=0 

( 



V 



17sin(ril) „2A:+2p 
r-C/n 



{2k + 2i)\tll{e) {2k + 2iy.tll-^\9)Ef 



+ 







k-\-l k-\-l+l 

um'-&] n [(2g-i)2-(e)2] 

q=0 q=l 



(2.45) 



We observe from this expression that type II defects will filter out all multiples of uj. For 
the pulse being once again of integer cycle length, this reduces to 



oo p 



Tn{2nu;,e,U„Eo) = J] E(-l)'"'"|fS (^^)'^' ' 

k,p=0 1=0 \ , 



2k + 21 
k + 1 — n 



(2.46) 



and 



oo p 



Tn{{2n-l)u^,e,Eo/e) = ^ ^^(-i; 



l+n+l *2fc 



2p+l 



{0) 



k,p=0 1=0 



22l-2p+l 



\k+p 



2p+l 



p\ {2k + 2l)\{2n-l)E^ 
I I {l + k-n+iy.{l + n + ky. 



(2.47) 



which are again independent of r. We observe that in this case we can not combine the Eq 
and UJ into a C/„. 



c 

CD 

O 
O 



2.5.3 One particle approximation 

In spite of the fact that we are dealing with a quantum field theory, it is known that a one 
particle approximation to the Dirac equation is very useful and physically sensible when 
the external forces vary only slowly on a scale of a few Compton wavelengths, see e.g. [43|. 
We may therefore define the spinor wavefunctions 



*j>,e(a;,t) : =ijf{x,t) 



*j-„,e(a;,t)t : = V/(x,t)^ 



a]{e) 



27r2p0 



2^2p0 



'tpj -x 



27rp0 



-ipj -x 



-.u,{e) 



27rp] 



(2.48) 
(2.49) 
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With the help of these functions we obtain then for the defect system 



(2.50) 



and the same function with u ^ v. Since this expression resembles a free wave, it can 
not be normalized properly and we have to localize the wave in form of a wave packet 
by multiplying with an additional function, g{p,t) in ( p. 20 ) and its counterpart g{x,t) in 
( |2.50 ), typically a GauBian. Then for the function <^f^g{x,t) = g{x,t)^f^ g{x,t), we can 
achieve that ||$|| = 1. 

2.5.4 Harmonic spectra 

We are now in the position to determine the emission spectrum for which we need to 
compute the absolute value of the Fourier transform of the dipole moment 



dt (*j>,e(x, t)tx*^^„_e(x, t)) exp int 



(2.51) 



We localize now the wave packet in a region much smaller than the classical estimate for 
the maximal amplitude the electron will acquire when following the laser field. We achieve 
this with a GauBian g{x,t) = exp(— x^/A), where A <C eEQ/uj'^. 

2.5.5 An example: Impurity of energy operator type 

As mentioned this type of defect, i.e. T){il),il:i) = gil)'ilj[x) can be obtained in a limit from 
the defect discussed in section 2.3.1. Coupling the vector potential minimally to it yields 



(2.52) 



by invoking the equation of motion. We can now determine the reflection and transmission 
amplitudes as outlined above 



Ri{9,g,A/e,y) = Ri{6,g,-A/e,-y) = R,{e, g, A/e, -y) = R,{9, g, -A/e,y) 

[yi - cosh 0] 6-2*2''^*°'^^ 
[1 - yA cosh 6*] - if + 1 + ^2 _ y2^2] sinh 9 ' 



(2.53) 



We denoted the differentiation with respect to time by a dot. The transmission amplitudes 
turn out to be 



Ti{9, g, A/e, y) = fi{9, g, -A/e, -y) = T,{9, g, -A/e, y) = f,{9, g, A/e, -y) 



1 _ y2^2 ^ _ |)2 



sinh( 



-[l-yA cosh 9]-i[^ + 1 + - y'^A'^] sinh 9 



(2.54) 



Locating the defect at ?/ = 0, the derivative of A does not appear anymore explicitly in 
( |2l53|) and (|2l5^ ), such that it is clear that this defect is of type I and admits an expansion of 



c 

CD 

O 
O 
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the form ( 2.41| ). With the explicit expressions ( ^l53|) and (|]5|) at hand, we can determine 
all the coefficients t2k{G) in ( |2.41| ) analytically. For this purpose let us first bring the 
transmission amplitude into the more symmetric form 



\T,i9,g,A/e) 



aoiO,9)+a2iO,g)A^ + ai{e,g)A'' 
ao{9,g)+a2{9,g)A^ + a^ie,g)A^-' 



with 



ao{9,g) = I6g^ + {4 + g^)^smh^9, 
a2{9,g) = 2g\4 + g^)smh^9, 



ao{e,g) = {g^ -4fsmh^9, 
o-i{(^-,g) = S'^sinh^^. 



(2.55) 



(2.56) 
(2.57) 



We can now expand \T{9^g^A)\ in powers of the field A{t) and identify the coefficients 



t2k{G,g) in ( 2.41 ) thereafter. To achieve this we simply have to carry out the series expan- 
sion of the denominator in ( |2.55 ). The latter admits the following compact form 



Y,<^2k{e,g)A^\ 



k=0 



aoie,g) + a2{9,g)A^ + ai{9,g)A^ 
with co{9,g) = l/ao{9,g) and 

C2k~2{0, g)a2{0, g) + C2k-i{0, g)a4,{9, g) 



(2.58) 



C2k{0,g) 



ao{0,g) 



(2.59) 



for > 0. We understand here that all coefficients C2k with /c < are vanishing, such 
that from this formula all the coefficients C2k may be computed recursively. Hence, by 



comparing with the series expansion (2.41), we find the following closed formula for the 
coefficients t2kiS,g) 



t2k{0,g) = [aQ{9,g) - ao{9,g)]c2ki9,g) k > 0. 
The first coefficients then simply read 

ao{9,g) 



to(.0,g) 



ao{0,g) 



mg) = ^-^\l 



\T{9,Eo = 
-to{0,g)] 



U{0,g) 



ao{6,g) 
ai{0,g) _ a2{0,g) 
a2{d,g) ao{9,g) 



8g'^i4 + g^)smh^29 
~ (16g2 + (4 + c/2)2sinh2^)2' 

t2{0,g), 



(2.60) 

(2.61) 
(2.62) 
(2.63) 



and so on. It is now clear how to obtain also the higher terms analytically, but since they 
are rather cumbersome we do not report them here. 



Having computed the coefficients t2k, we can evaluate the series ( 2.42 ) and ( 2.45 ) in 
principle to any desired order. For some concrete values of the laser and defect parameters 
the results of our evaluations are depicted in figure 3. 

The main observation from part (a) is that the defect acts as a filter selecting higher 
harmonics of even order of the laser frequency. Furthermore, from the zoom of the peak 



c 

CD 

O 
O 
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Figure 3: Fourier transform of the transmission probability for a single (a) and double (b) defect 
with Eq — 2.0, g = 3.5, 9 = 1.2, oj = 0.2. Harmonic emission spectrum for a single (c) and double 
(d) defect with Eq ^ 2.0, g = 3.5, 9 = 1.2, ui = 0.2, A 6. 



regions, we see that there are satellite peaks appearing near the main harmonics. They 
reduce their intensity when r is increased, such that with longer pulse length the harmonics 
become more and more pronounced. We also investigated that for different frequencies cu 
the general structure will not change. Increasing the field amplitude £"0, simply lifts up the 
whole plot without altering very much its overall structure. We support these findings in 
two alternative ways, either by computing directly ( 2.39| ) numerically or, more instructively, 
by evaluating the sums (|2l^ ) and (|2l^) . 

Part (b) shows the analysis for a double defect system with one defect situated at 
X = and the other at x = y. The double defect amplitudes are computed directly from 
(2.12) and ( |2.13| ) with the expression for the single defect ( p. 53 ) and ( |2.54 ). Since now 
both A and A appear explicitly in the formulae for the i?'s and T's, it is clear that the 
expansion of the double defect can not be of type I, but it turns out to be of type II, i.e. 
of the form ( |2.44| ). Hence, we will now expect that besides the even also the odd multiples 
of uj will be filtered out, which is indeed visible in part (b) for various distances. Here 
we have only plotted a continuous spectrum for y = 0.5, whereas for reasons of clarity, 
we only drew the enveloping function which connects the maxima of the harmonics for 
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the remaining distances. We observe that now not only odd multiples of the frequency 
emerge in addition to the ones in (a) as harmonics, but also that we obtain much higher 
harmonics and the cut-off is shifted further to the ultraviolet. Furthermore, we observe a 
regular pattern in the enveloping function, which appears to be independent of y. Similar 
patterns were observed before in the literature, as for instance in the context of atomic 
physics described by a Klein-Gordon formalism (see figure 2 in p4[). 

Coming now to the main point of our analysis we would like to see how this structure 
is reflected in the harmonic spectra. The result of the evaluation of ( ^.51 ) is depicted in 



figure 3 parts (c) and (d). We observe a very similar spectrum as we have already computed 
for the Fourier transform of the transmission amplitude, which is not entirely surprising 
with regard to the expression ( |2.51| ). The cut-off frequencies are essentially identical. From 
the comparison between X and the enveloping function for T we deduce, that the term 
involving the transmission amplitude clearly dominates the spectrum. 

The important general deduction from these computations is of course that harmon- 
ics of higher order do emerge in the emission spectrum of impurity systems, such that 
harmonics can be generated from solid state devices. 

3. Conductance from the Kubo formula 

Having characterized various features of defects, I will proceed with the main theme of the 
talk, that is the computation of the DC conductance. In the absence of impurities it can 
be obtained from the Kubo formula in the form 

1 f°° 

G(r) = -lmi^ J Jte^-' {J{t)JmT,m- (3.1) 



We proposed in |jl| a generalization of ( p.l| ) in the form of (1.1). The key quantity needed for 
the explicit computation of (^) or (|l.l| ) are the occurrence of the temperature dependent 
current-current correlation functions (J(r) J(0))j.^ or ( J(r) J(0))y ^, respectively. 

In the zero temperature regime two-point correlation functions can be computed in 
general by means of the form factor bootstrap approach ||l5|, |l^ . In this approach one 
expands the two-point function between two local operators O and O' in terms of the series 



g -rm-i cosh 9i 



n=lfii---fi„ ^ ' i=l 

xFf l^-^"(0i, . . . , e„) [i^f >-'^"(ei, . . . , On)] * , (3.2) 

where we choose = (— ir, 0). The form factors are defined as matrix elements of the 
local operator 0{x) located at the origin between a multiparticle in-state and the vacuum, 

Ffi^-^"(^i,^2...,^n) := {o\omzl{ei)zl{e2)...zl{en)). (3.3) 



The expansion (3^) is simply obtained by inserting complete states on the r.h.s. One may 



proceed similarly by inserting one more set of complete states when a defect is present and 
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obtains 

(J(r)Z„J(0))^^o = 



f dOi--- dOndOl ■■■dOm c.J|;,i...m„ (g n ^ 

-2^ 2^ y m!n!(27r)"+'» " {^i---»n} 

n,m,=l tl\---^ln\l'l-- Um 



x(z^„(0„)...Z^,(0i)|Z,|Z,,(^i)...Z,^(e„))F^I— -(ei...0„)*e -1 . (3.4) 

This means there are three principle steps left in order to obtain the conductance from 
the expression in (|1 . 1| ) . (a) The computation of the form factors (|3.3| ) and the matrix 
elements involving the defect operator occurring in ( |3.4D . (b) The integration in r and 
(c) the limit lo ^ 0. Step (a) can be performed in two alternative ways either by solving 
certain consistency equations for the form factors and defect matrix elements or by direct 
computation. For the latter we require a representation for the particle creation operators 
Z^{6), the defect operator Zq, and the local operator 0{r) which is the current in this case. 

3.1 The massless limit 

Remarkably when carrying out the massless limit of the above expressions, the steps (b) 
and (c) can be carried out generically. To perform such a limit we proceed according to 
the massless limit prescription as suggested originally in ^5|. It consists of carrying out 
the limit m — > in the high energy regime. In order to do this one replaces in every 
rapidity dependent expression 9 hy ziz a, where an additional auxiliary parameter a has 
been introduced. Thereafter one takes the limit a oo, m ^ while keeping the quantity 
m = m/2exp(o") finite. For instance, carrying out this prescription for the momentum 
yields p± = ibmexp(ib0), such that one may view the model as splitted into its two chiral 
sectors and one can speak naturally of left (L) and right (R) movers. For the form factors 
in (|3.4| ) the massless limit yields 

hm Ff + lJ^a, ...,6^+ rj^a) = • • • , ^n), (3.5) 

with r]i = ±1 and Ui = R for rji = + and Vi = L for rji = —. Namely, in the massless limit 
every massive n-particle form factor is mapped into 2" massless form factors. Using these 
expressions, performing a Wick rotation and introducing the variable E = XliLi ^j^^S 'we 
obtain from 1^2 



oo 



n,m=l /ii---Atn;i^i--i^m ^ ' 



z^jOn) . . . z^^{e^)\z^\z^^{e,) . . . z^ie^)) p-^^^-^^ih. ■ ■ . ,e~™)*e--^. (3.6) 

We note that for the massless prescription to work, the matrix element involving the defect 
Za can only depend on the rapidity differences, which will indeed be the case as we see 
below. Performing the variable transformation 9^ — > In E'/rhn — XlILi '^i/™'"^^*' 
write the r.h.s. of (13.61) as 



E 

n,m=l ^\---^n\^'\---Um 



J n!(27r)" J m!(27r)' 

-oo — oo 
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-irE' 



(3.7) 



We substitute now this correlation function into the Kubo formula, shift all rapidities as 
di ^^E'/rhn, Oi 6i+ In E'/rhn, use the Lorentz invariance of the form factors^ 

and carry out the integration in dE' 



lim 



2s-2 



/ii---/i„;i/i---i/„ 

n-1 



dOi--- dO. 



71-1 



dOi--- d9„ 



1 



n!(27r)" J m!(27r)- i - ^^-^ rm/mne^^ 
Z^JHl - Y.Zi m/mne'^)) . ■ ■ Z^^{ei)\Z^\Z^^{h) ■ • • Z^JO,,)) (3. 



n-l 



,ln(l-^m,/m„e^O)F^|' 



J\ui...Un 
■Rm 



^y 



i=l 



We state various observations: Since the matrix element involving the defect only depends 
on the rapidity difference, it is not affected by the shifts. Operators with Lorentz spin s = 1 



play a very special role in (3^), which makes the current operator especially distinguished. 
In that case the r.h.s. of Q becomes independent of the frequency uj and the limit is 
carried out trivially. Furthermore, since the final expression has to be independent of m„, 
we deduce that the form factors have to be linearly dependent on m„. 

3.2 Realization of the defect operator 

A realization of Z^ can be achieved very much in analogy to a realization of local operators, 
i.e. as exponentials of bilinears in Zamolodchikov-Faddeev operators For the case of 
a boundary a generic model independent realization for the boundary operator B was 
originally proposed in for the parity invariant case, i.e. R = R . This proposal was 



generalized to the defect operator in [26| with the same restriction and for self-conjugated 
particles. Here we extend this realization in order to incorporate the possibility of parity 
breaking as well as non self-conjugated particles. A non-trivial consistency check for the 
validity of our proposal will be ultimately provided when exploiting it in the computation 
of the conductance, obtained by entirely different means as will be presented in part II. 
The realization we want to propose here is a direct generalization of the one presented in 
p^ , namely 

Z„=:exp[— / Da{e)de] : , (3.9) 
where : : denotes normal ordering and the operator Da{9) has the form 



(3.10) 



^Denoting by s the Lorentz spin of the operator O and A being a constant, the form factors satisfy 



CD 

o 
o 



F' 



0|a«i...m„ 



(Si + A, . . . , 0„ + A) = e"- F, 



S\ rpO\l-Ll...fln 



{Oi,. 
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with K-{d) := RfCf - d), kf{e) := i?f (f - 9), W^O) ■= T^'f - 0) and := 
T"{^ — 0). In comparison with p6| we have used a sUghtly different normalization factor, 
since in general we have contributions in the sum over i in ( 3.1C| ) including both particles 



and anti-particles, as for the complex free Fermion we shall treat below. Following the 
arguments given in |25], the operator Da{0) depends on the amplitudes R{0), T(6), R{0) 



and T{9) with their arguments shifted, as considered also in |24 



3.3 Defect matrix elements 

Having now a concrete generic realization of the defect ( |3.9| ), we can compute the defect 
matrix elements. One way of doing this is to solve a set of consistency equations which 
relate the lower particle matrix elements to higher particle ones, similar as in the standard 
form factor program [15, 16, 17]. Such kind of iterative equations were proposed in |^] for 



a parity invariant defect and for a real free fermionic and bosonic theory. We generalize 
this here and note first that the operator (p. 91) becomes 



1 

lim Z„=:exp[— / dOj^ Zj{9)Z,{e)] (3.11) 

and the defect should act in this case as the identity operator, which fixes our normalization 
to {Zi[9\)ZaZ^-(9'2)) = 2t: 5{9i2)5ij after having contracted according to Wick's theorem. 



3 

For two particles we find 



{Z,{9i)Zi{92)Z^) = vrETf (02)5(^12), (3.12) 
{Z^ZI{9^)ZI{92)) = 7TK^{9ir6{9i2), (3.13) 
{Z,{9i)Z^z}{92)) = 7rW,"{9i)6{9u)Sij. (3.14) 

For later convenience we have introduced the functions 

kr{9) = Kf{9) + Sn{-29)K-^i-9) = kf{9) + Ss{29)k^{-9), (3.15) 
Wt{9) = W^{9) + W^{-9y = W-^{-9) + W-^{9y = W-^{9)\ (3.16) 

since the Kf, isr",Wj° and amplitudes defined before will repeatedly appear in the 
combinations ( p. 15]) , ( |3.16| ) in what follows. The latter equalities in ( 3.15| ), (|3.16| ) follow 
simply from 

Wr{9) = Wr{-9) = W-^{iTT-9)\ kf{9) = S^,{29)Kf{-9) = Sij{29)k^{iTT-9)\ (3.17) 

which are in turn consequences of the crossing-hermiticity properties ( |2.3| )-( |23^ . With 
these matrix elements we can construct the ones involving more particles recursively from 

,,,e,^e[... 9'^) := [z,M ■ ■ ■ ^Mi(^i)^- zli^'i) ■ • • zIM) = 

n ■••'''■••''''^"■""(^m ...9i...92,9[...9'J (3.18) 

1=2 p=l 

n l-l 

- n •■■'^^''^■••''■•""(^™ ...92,9[...9l... 9'J 
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i-i 



n 

1=2 p=l 

m l—l 

1=1 p=i 



(3.19) 



fll...fJ.lU2..M„ I 



n - ■ ■ 



90- 



Here we denoted with the check on the rapidities 6 the absence of the corresponding particle 
in the matrix element. It is clear from the expressions ( |3.9D and ( |3.10| ) that the only possible 
non- vanishing matrix elements ( |3.18| ) are those when n + m is even. Taking ( |3.12 )-( 3.14| ) 
as the initial conditions for the recursive equations ( |3.18| )-( |3.19| ), we can now either solve 
them iteratively or use (^]^) and evaluate the matrix elements directly. Closed solutions 
for these equations have been presented for the first time in Q. 

3.4 Free Fermion wire with impurities 

At this point we have to abandon the general discussion and consider a concrete theory, 
which for the reasons already explained we choose to be the complex free Fermion. Then 
the generators of the ZF-algebra Zi(9), z\{9) are just the usual creation and annihilation 
operators ai{9), a\{9). 

3.4.1 Defect matrix elements 



Let us now use (p.9|)-(3.10) in order to evaluate matrix elements involving the defect oper- 
ator. In what follows, the most relevant matrix elements are those involving four particles, 
for which we compute 

{a0i)a,{92)Z^a\{9:,)a\{9^)) = wf,{9i92)S{9u)S{923) + kfi{9i9^)d{9i2)d{93A), 
{ai{9i) ai{92)Z^ a]{9:,) a]{9^)) = -tt^W^{9i)W^{92)5{9M92a%j, 



{ai{9i)ak{92)a,{93)Z^a\{9^)) = 7:''Wt{9i)K^{-92) [6{9u)6{923) - 6{9i2n9si) 

{a,{9i)Za4{92)al{93)4{94)) = 7T^Wt{9i)k^{-93r [s{923)S{9u) - 6{9i2)d{9u) 

with the abbreviations 

wU0i92) = 7r'Wt{9,)Wr{92) and A;-(0i,e2) = ^'i^f (^i)i^r(^2)* • (3.20) 

One can now try to find solutions for all n-particle form factors either from ( p.l8| )-( |3.19| ) or 
by direct computation. For instance for the stated choice of particles involved, we compute 



iki 



77m x(h) nx(ii) / n 



'2n) 



min(n.m) 

E 



(-1) 



m+n—2k ^n+m 



k=0 



{m - k)\{n - k)\k\k\ J_ 



X det A^'^iPl ...(32n;9[... 9'2J det ^''"(/^Sn+l • • • (32n+2m; 9l . . . 92m) 



X J] W^{(52p)W-^{l32p-lW2p - P2n+2p)S{(32p-l " /?2n+2p-l) 
p=l 

n n+m 

X II Kr(lh.r,Vd[lh.r, + li-?.v-i) n Kf{l32p)5{(32p + f52p-i) 



d(3i . . . d(52n+2m 



(3.21) 



n Kf{(i2prmp+^ 

p=l+k 



c 

CD 

O 
O 



p=l+k+n 
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where A^{9i . . .9f,9[ . . . 6'^) is a rank £ matrix whose entries are given by 

Ai = cos^\{i-j)7:/2]6{ei-e',), l<i,j<i. (3.22) 

and references therein. Likewise we 



.,^-cos'[{i-j)n/2]6{e,-9'^) 

The matrix elements are computed similarly as in 
compute 

7r"(-l)"-i 



rpnxi+mxi I 
a V 



C„ . . . C7l, (7^ 

xdetS" 



9' ) 

'ml 



d[3i...d(3nj[wt{ek) 



k=l 



„ . . . (71, pi 



;/3i.../3„)det^"(/?i.../3„;0'i...O, 



(3.23) 



where we introduced a new rank I matrix B^{9i . . . 6f,9'i . . . 0'^) whose entries are now simply 
given by 

Bl^ = 5{ei-e]), i<i,j<e. (3.24) 

One can verify explicitly ||^] that these expressions indeed satisfy (3.18) and ( 3.19| ) . 



3.4.2 Conductance in the T = m = regime 

It is well-known that for a free Fermion theory (also for a single complex free Fermion) the 
conformal C/(l)-current-current correlation function is simply 



(J(r)J(O)) 



1 



T=m=0 



(3.25) 



This expression can also be obtained by using the expansion (3.2), together with the mass- 
less prescription as outlined above and the expressions for the only non-vanishing form 
factors of the current operator in the complex free Fermion theory 



pj\ti 



(9,9) 



-F/I"(e, 9) 



e+0 

-iirme 2 



(3.26) 



± 2 — --t 2 

In particular, the massless limit of the previous expressions gives, according to the massless 
prescription. 



Fi^S{9,9) 
Fit {0.0) 



F'l^S^OrO) 



8+e 
-27:1 me 2 



(3.27) 



Fi\^{9, 9) = F^^l\9, 9) = Fif{9, 9) = fI\^\9, 9) = F'^f{9, 9) = Q. (3.28) 

We these expressions we can evaluate (U) to (|^). We may the insert ( |3.25D into ( |3.1] ) 
and the problem is reduced to find the Fourier transform of the function r~^, which is given 
by f^oo e*'^'"r~^ = —TTLJ for u) > 0, with V denoting the principle value. This yields in 
the absence of a defect G{0) = 1/2tt, in complete agreement with the well-known classical 
expression for the conductance in a wire without any impurities, see for instance (47|. 

For the more complicated situation of n defects Za^ ■ ■ ■ Za„ located in space at positions 
Uai ■ ■ ■ Van ) compute in the zero temperature and zero mass regime 



(J(r)Z,,-..Z„„J(0)) 



~ 2 
m 



T=m=0 



+ 



E 



d9, 



'2rm cosh 81 



kf''{9,) 



d9i 



9\ —rme"i 



w: 



a\R, 



7^) / '^^e2-rrhe02^a\R 



{92 



(3.29) 



c 

CD 

O 
O 
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The functions W°'^^{6), K"^^{6), . . . defined in ( p.29D are the massless hmits of the cor- 
responding functions W°'{9), Kf{6), ... For all the defects we considered, it turned out 
that the first contribution to the previous correlation function is actually vanishing, so 
that (3.29) is considerably simplified. In many of the examples, this is due to the fact 



that the amplitudes Kf{0) are vanishing in the first place, as a consequence of the cross- 
ing relations ( ^.17 ). The vanishing of the reflection part in ( 3.29| ) also occurs in some 



cases as a consequence of the parity of the function Kf{6). For instance, we find that, 
for the energy operator defect such function, although initially non-vanishing, satisfies 

kf{e) = -K^i-e), such that limm^oj^^dek^ie)* = o. 

We can now either use ( |3.29|) to compute the conductance or evaluate the expression (^^) 
directly in which the frequency limit is already taken, in both cases we obtain 



= 2(^^ / dO e' wf^'^'^Ml - e')e] . (3.30) 

* -oo 

There are, in addition, further generic results which can be obtained independently 
of the specific form of the defect. We present them at this stage and will confirm their 
validity below by some specific examples. Specializing to the case in which all i defects are 
of the same type and equidistantly separated, i.e. y = yai = ■ ■ ■ = ya„- We can identify 
two distinct regimes 



<l««(0i e,) = vr^ I wfy.Wf^'ie.Y for finite y 

for y^O ^ ' 

where we used in addition ( ^.161) . Supported by our explicit examples below, we find that 



for y ^ in (l3.3l[) the amplitudes W"^^{9) become independent functions of the rapidity. 



As we have already argued above 



kT''{Gi,Q2) = 0. (3.32) 



It will turn out, that the two regimes specified in (3.31) are also of a very distinct nature 
in the TBA context as presented in part II. 

3.4.3 A wire with impurities of energy operator type 

Let us exemplify the working of the above formulae with a concrete defect operator. As a 
simple example we choose the energy operator defect as presented in section 2.3.1. Con- 
sidering first a wire possessing a single defect of this type, we compute 

,t,^.„s 4 cos B cosh^ 9 . 2i sinSsinh^ a\RR,„ „ „^2 

Wi{9) = — — — — , K^[9) = - — — , w,j [9i92) = [2iTcosBy 

* ^ ^ cosh 26* + cos 2S * ^ ^ sm B - cosh 9 v . / v ^ 

(3.33) 

with B being the effective coupling constant as defined in the caption of figure 1, such that 



(J(r)Z„ J(0))^_^o = ^ =^ G"(0) = (3.34) 
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It will turn out that this is in complete agreement with the corresponding result from the 
Landauer formula ( |1.1| ). 

Proceeding in the same way for a wire with two or four impurities we evaluated [P in 
the regime y ^ r 



=m=0 

{J {r)Za^Za2ZoizZaiJ {^)) T=m=a ' 

In the regime y ^ 0, we obtained [ffl] 



4 [1 + sin^ B] 
r2 [cos2(25) -3] 
2 1 + sin^ B 
vr[3-cos2(25)] 



2 ' 



2 ' 



1 

1 



1 + 



1 + 



cos'^B 



[cos^B- 2(1 + sin2s)2]2 
cos^5 



47r V [cos4 B-2{1 + sin"^ B)2]2 



lim (J {r)Z a, Za2 J {0))T=m=0 

lim G"i"2(o) 



hni (J (r ) Z^.^ J (0) )-r^^^o 

y^O 



1 COS^ 5 

^(l + sin2 5)2 ' 
1 cos^ B 
2^(1 + sin2 5)2' 

1 / cos'^S 



cos^B- 2(1 + sin2s)2 
cos^ i? 



lim G°i°2a3a4(o) 



1 



(3.35) 
(3.36) 
(3.37) 
(3.38) 

(3.39) 
(3.40) 
(3.41) 
(3.42) 



2tt Vcos^S - 2(1 + sin2 5)2. 

It will turn out that we can reproduce these expressions by evaluating the Landauer formula 
( |1.2| ) when computing the densities with the help of the TBA. This will now be outlined 
in part II together with the general conclusions concerning also this part. 
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